Light Hadron Spectrum in Quenched Lattice QCD with Staggered Quarks 
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Without chiral extrapolation, we achieved a realistic nucleon to p- meson mass ratio of mjv/m p = 
1.23 ± 0.04(statistical) ± 0.02(systematic) in our quenched lattice QCD numerical calculation with 
staggered quarks. The systematic error is mostly from finite-volume effect and the finite-spacing 
effect is negligible. The flavor symmetry breaking in the pion and p meson is no longer visible. 
The lattice cutoff is set at 3.63 ± 0.06 GeV, the spatial lattice volume is (2.59 ± 0.05 fm) 3 , and 
bare quarks mass as low as 4.5 MeV are used. Possible quenched chiral effects in hadron mass are 
discussed. 



Reproducing the known light hadron mass spectrum 
is the most important test the numerical lattice QCD 
is yet to pass, in spite of the steady progress [] since 
the pioneering works by Weingarten and Hamber and 
Parisi []. The main obstacle is the difficulty in including 
light dynamical quarks, and consequently the available 
full-QCD calculations still suffer from too heavy quark 
mass, too coarse lattice spacing or too small lattice vol- 
ume []. On the other hand, with the quenched approx- 
imation where one neglects dynamical quark loops, re- 
cent calculations use small enough lattice spacing and 
large enough lattice volume to understand the system- 
atic errors arising from them []. Indeed recent quenched 
calculations [] collectively have shown that both of these 
errors are smaller than the statistical noise, albeit with 
rather heavy quarks. Yet these calculations left three 
major problems: nucleon to p- meson mass ratio is too 
high, pion to p-meson mass ratio is too high, and ex- 
trapolating the results to more realistically light quark 
mass values is necessary but difficult because of the sub- 
tle issue of the quenched chiral effect [] . Hence quenched 
calculations with realistically light quark mass values on 
a large enough and fine enough lattice are desirable. 

Helped by the results from ref. [], we choose a gauge 
coupling of 6/g 2 — 6.5 and a lattice volume of 48 3 x 64. 
We will find later in this letter that these parameters 
correspond to the lattice spacing of a = 0.0544(9)fm or 
the cutoff of a -1 = 3.63(6)GeV and a spatial volume of 
(2.59(5)fm) 3 . We use staggered quarks because it is defi- 
nitely superior to the Wilson one in controlling the quark 
mass and hence in investigating the issue of quenched chi- 
ral effect: the quark mass is well defined and protected by 
the remnant U(l) chiral symmetry in the former while in 
the latter one encounters a difficult problem of defining 
the critical hopping parameter under the inevitable pres- 
ence of exceptional gauge configurations [] . For the gauge 
part we use the single-plaquette Wilson action because 
our lattice spacing is fine enough. 

A combination of multi-hit Metropolis and over- 
relaxation algorithms is used to generate Monte Carlo 
samples of quenched gauge field. Separation be- 



tween propagator sampling is 2000 such updates: 1000 
Metropolis interleaved with 1000 over-relaxation. This 
is proven good enough from the auto-correlation analy- 
sis of the obtained pion propagators. We use the con- 
jugate gradient (CG) method for inverting the stag- 
gered quark Dirac matrices. A few different sizes of 
corner and even wall sources with bare mass values of 
m q a = (0.05,0.04,0.03,0.02,0.015,0.0075) (Set I) and 
m q a = (0.01,0.005,0.0025,0.00125) (Set II) are used 
for calculating staggered quark propagators. These two 
propagator sets are obtained from two almost indepen- 
dent sets of gauge field configurations: they share only 
a few gauge field configurationsin common. This is to 
further reduce correlations and should lead us to better 
comparison of fitting results from one set with those from 
the other set. They are all combined with point sinks. 
Set I and Set II together, the bare quark mass varies 
for over a factor of 40 and provides us a good theater 
in studying the chiral behaviors. We tried two different 
kinds of wall definition and a few different source sizes 
to eliminate systematics arising from using a single kind 
and size. 118 of quark propagators is collected for Set 
I, and 250 for Set II. The numerical algorithms are ba- 
sically the same as in ref. [] and technicalities associated 
with our implementation on VPP-500 vector-parallel su- 
percomputer is given in ref. [] . Preliminary reports of the 
obtained results were given in ref. [] . 

Tabic H summarizes our estimates for the pion, p- meson 
and nucleon mass values. Here m^a is the mass estimate 
for the mass of Goldstone pion, while m^a is the es- 
timate for non-Goldstone pion extracted simultaneously 
with the estimate for its parity partner scalar /o/ao- Sim- 
ilarly m p a is the estimate from the vector meson part- 
nered with b\ axial while m P2 a is the estimate extracted 
simultaneously with its parity partner a\ axial. The nu- 
cleon mass to jv a is from the even-source results which 
gave better signals than the corner-source ones. Fit- 
ting is done by minimizing the correlated x 2 calculated 
from a single elimination jack-knife data set. From these 
data a few immediate conclusions follow: 1) In Fig. [j] 
we show our Edinburgh plot, m^r/m^ vs. m^/rrip. For 
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our lightest bare quark mass of m q a = 0.00125, we get 
m N /m p = 1.230 ± 0.035 and m^/m p = 0.273 ± 0.006. 
The former is in good agreement with the observed value, 
although the latter is still about 50% larger. 2) The fla- 
vor symmetry breaking estimated by (m^a — m,«) and 
(m P2 a — m p a) is generally small. It decreases as we de- 
crease the bare quark mass so much as to be eventually 
hidden below statistical errors. The symmetry is restored 
well enough. 

Now let us turn our attention to what we can learn 
about the systematic errors. First, the finite volume ef- 
fect: from Table we see that m p a between m q a — 0.0025 
and m q a — 0.00125 does not change within statistical 
error. Assuming, rather safely, that the p-meson mass 
dependence on the quark mass is mild in this region 
around m q a — 0.00125 and taking the calculated re- 
sult as physical, we estimate the physical size of our lat- 
tice cutoff a -1 be m p (physical) / m p (m q a = 0.00125) = 
0.7700(8)/0.212(4) « 3.63(6)GeV~ 1 or a = 0.0544(9)fm. 
It follows that our lattice has spatial extent of about 
2.59(5) fm. Using nucleon mass instead of p-meson mass 
results in a consistent estimate. It has been argued that 
the lattice-QCD finite-volume effect is sensitive to the 
pion Compton wavelength to" 1 on the lattice: when it is 
large compared with the lattice, the finite volume effect 
is expected to fall like 1/V [], and like exp(— m^L) other- 
wise [] . An extensive study made by the MILC collabora- 
tion [] report that the finite- volume effect on a (2.7fm) 3 
lattice is smaller than 1 %. This should translate into at 
most 1.3 % effect for our lattice volume of (2.59(5)fm) 3 
when 1/V dependence is assumed. With the exp(— m^L) 
dependence the effect is smaller than this except for the 
case of our lightest bare quark mass of m q a — 0.00125, 
where the effect is expected to be slightly larger with 
TOjrL = 2.76(4). We also compare our data with existing 
data at the same 6/g 2 = 6.5 but on a smaller 32 3 x 64 
lattice [] for m q a = 0.01, 0.005 and 0.0025. Here we see a 
0.7±1.0 % (m q a = 0.01) to 4.4±3.2 % (m q a = 0.0025) ef- 
fect in TOpCt (heavier m p a on a larger lattice), and 4.7±1.0 
% (m q a = 0.01) to 6.3 ± 3.2 % (m q a = 0.0025) effect in 
TOjva (lighter m^a on larger lattice). These are con- 
sistent with 4.4 % effect expected for the 32 3 volume 
assuming the 1.3 % effect on the 48 3 volume and the 
1/V behavior for 32 3 volume. Therefore we estimate the 
finite-volume effect in the current nucleon to p-meson 
mass ratio /m p result is V2 x 1.3% x 1.23 ~ 2.3%. 

Finite-spacing effect: here we expect 0{a 2 ) flavor 
breaking effect among various definitions of staggered pi- 
ons and p mesons. However as we already discussed, the 
breaking is hardly visible in our data alone. Comparison 
with earlier works [] at lower values of Q/g 2 (< 6.2) re- 
inforces this observation. Finite lattice spacing effect in 
the mass ratio like tojv/to p should be even smaller than 
that in m p and tojv individually. Though there are po- 
tential O(a) effects to TOat /m p from the flavor symmetry 
breaking in m^a [], it should be negligible as the breaking 
in m^a is already hardly visible. 



In addition to these systematic errors mentioned in 
the above, we considered whether the size of the quark 
field wall source introduces a systematic bias in choosing 
the best fit for hadron mass (depending on the size of 
wall source, excited hadronic states can couple to the 
wall differently []). For the three different wall sizes, 
12 3 ,24 3 , and 32 3 , we gathered 300 hadron propagators 
with m q a = 0.01 and 124 hadron propagators with 
m q a = 0.00125. Fig. |2| shows the effective mass for pion 
from these different source sizes at m q a = 0.01. We see 
an expected tendency that the plateau sets in later for 
smaller sources. More importantly, within the statistical 
error they eventually agree with each other and elimi- 
nates ambiguity in defining a plateau. Other effective 
mass plots exhibit the same behavior and help defining 
plateaus. Thus, our choice of the best fit is less biased 
by a subjective choice of a plateau. 

Thus we established a good enough control of the 
systematic errors arising from the finite volume, finite 
cutoff, and choice of plateau in the effective mass, to 
start discussing the quenched chiral log problem. In 
Fig. ^, m 2 /m q vs. m q a is plotted. The fitting form, 
m 2 , — Cq + C\vn q + C2IT1 2 , has been tried on Set I and 
Set II separately with correlations among different to^'s 
included. This form has been suggested by the finite vol- 
ume effect on the pion mass []. At best we get Co = 
(1.9 ± 0.2) x 10~ 3 , d = 2.20 ± 0.02, C 2 = 14.4 ± 0.2 with 
confidence level (C.L.) of 1.05 xl0~ 5 (x 2 /d.o.f. = 8.6) for 
Set I, and C = (1.7±1.1) x 10~ 4 ,Ci = 2.51 ±0.04, C 2 = 
-(0.30 ± 0.33) with C.L. of 2.2 x 10~ 4 (x 2 /d.o.f . = 13.6) 
for Set II. Neither fit is satisfactory. In particular, the 
fit to Set I overshoots in the quark mass region cov- 
ered by Set II. This suggests our data in the small mass 
region are much less singular than the l/m q behavior 
of the finite volume effect. For Set II, a fitting form 
inspired by the quenched chiral perturbation theory [], 
log to^./ Tn q — c — <51ogTO 2 2 has been tried, where the 
correlations among to^ and m- n2 have been fully taken 
into consideration. Fitting all the data from Set II to- 
gether gives c = 0.815 ± 0.037 and 5 = 0.027 ± 0.010 
with C.L. of 9.7 x 10~ 4 (x 2 /d.o.f. = 6.93). This confi- 
dence level or x 2 /d.o.f. is better than that for the fitting 
form considered in the above but it is still marginal. Fit- 
ting m q a = 0.005,0.0025 and 0.00125 data only gives 
c = 0.691 ± 0.052, S = 0.057 ± 0.013 with C.L. of 
5.8 x 10~ 2 (x 2 /d.o.f. = 3.59) while fitting to the form, 
rn 2 = Co + CiTO g , for these three quark masses gives 
Co = (2.29 ± 0.98) x 10~ 4 , d = 2.47 ± 0.02 with C.L. 
of 1.6 x 10~ 4 (x 2 /d.o.f. = 14.2). One can try to add one 
more parameter to the above quenched chiral log fit, in 
order to fit all the data from Set II. Indeed, such a fit 
gives an improved C.L. (7.0 x 10~ 3 ) but similar modifi- 
cation (from Co + CiTO 9 to Co + C\m q + C2TO 2 ,) does not 
improve C.L. of the finite volume fitting form. Thus we 
think that the quenched chiral fitting form describes our 
data better and the quenched chiral logarithm behavior 
discussed in the to 2 / m q in [] is not a finite lattice volume 
artifact. 
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For nucleon and p- meson mass values, various fitting 
forms are suggested and tried in ref. [], where origins 
of each terms in these fitting forms were discussed. Al- 
though these forms should be considered only after the 
continuum limit is taken, we try these fitting forms on 
our data since finite lattice spacing effect on our data is 

1/2 

small as we discussed in the above. m q and m q log m q 
terms are from quenched chiral perturbation theory con- 
sideration, and m q (from tree level), m q ^ 2 , m 2 , m 2 \ogm q 
(from one loop correction) terms are present both in 
quenched and ordinary chiral perturbation theory. Fol- 
lowing them, we studied chiral extrapolation in m p a and 
77i a? a using hadron masses from Set I and compared di- 
rectly with those masses from Set II. Since Set II has 
four data points, we also tried fits to hadron masses from 
Set II when the number of fitting parameters is less than 
four. Correlation among hadron masses is included in 
the fitting. For Set I, among the twelve fitting forms sug- 

1 /2 

gested in ref. [] , we can definitely rule out a + bm q and 
a + bm q because the confidence levels are so poor. Fitting 

to a + brri^ 2 + cm q + drn'J 2 and a + bm q ^ 2 + cm q + dm 2 
return either d with error more or less equal to d sug- 
gesting that the d term is not necessary, or a nega- 

1/2 

tive d. On the other hand, the fit to a + bm q + cm q 
gives a positive value for b, which is inconsistent with 
quenched chiral perturbation theory []. Similarly, fits to 

a + bm q + cm^/ 2 + dm 2 , a + bm^ 2 + cm q + dm q logra 9 , 
and a + bm q + cm 2 + dm 2 log m q give either d with error 
as large as d or larger than d implying that fitting forms 
of a + &777,,j + cm q ^ 2 , a + brricj 2 + cm q , and a + bm q + cm 2 
are preferred, or a negative value for d. All the fits to 
a + bm q + cm q ^ 2 and a + bm q ^ 2 + cm q returns a negative 
c. The fit to a + bm q + cm q log 777. g gives a negative c. 
The confidence level of the fit to a + bm q + cm 2 log m q 
is 7.7 x 10~ 4 for m N a and 1.3 x 10~ 2 for m p a. None of 
these extrapolations agrees well with results from Set II. 
In contrast, all the fitting forms except for the a + bm q 
to Set II give high C.L., telling us that the associated 
statistical errors of Set II are too large. Following ref. [], 
we also tried fitting itin + Xj^m V2 and 77i p + \ p m 7r2 to 

a + bm q + cm q ^ 2 + dm 2 or a + bm q + cm' q 2 for various A's 
but none of them improves confidence level significantly. 

In Fig. ^, our m^a and 77i p a are plotted against m 7T2 . 
The leftmost circles show expected uiNd and m p a for 
physical quark mass which is deduced by the nucleon 
mass and the pion mass from experiment and by the lat- 
tice spacing, a, obtained in the above. Quenched chiral 
perturbation theory suggests a linear dependence on m l2 
with a negative coefficient. A fit gives a positive coeffi- 
cient with poor confidence level, 0(1O -4 ) although the fit 
looks good to eyes. Our data in Fig. || clearly show that 
the coefficient of linear term is small, in agreement with 
small 8 (0.02 ~ 0.06) from our m 2 /m q data. The small- 
ness of the coefficient for the linear m W2 term and the 
disagreement between the extrapolated results from Set 
I and the actual Set II data, lead us to conclude that one 



needs very high statistics to trust chiral extrapolation of 
hadron mass data when lattice calculation is performed 
with heavy quark masses in order to obtain a result sim- 
ilar to MILC collaboration []. It is interesting that for 
77i 9 a < 0.005, the obtained values of pion and p-meson 
mass would allow p — > tttt decay. Influence of such decay 
mode on p-meson mass needs further investigation in the 
context of quenched approximation. 

In conclusion our numerical calculation of quenched 
QCD with (2.59(5)fm) 3 spatial volume and 3.6 GeV cut- 
off yields a realistic mass ratio of m,^jm p = 1.230 ± 
0.035 ± 0.023 at the bare quark mass of 77i 9 a = 0.00125 
~ 4.5 MeV, where the first error is statistical and the 
second is finite-volume. Finite lattice spacing effect is 
negligible. This result is obtained without chiral ex- 
trapolation. Flatness of m/yci and m p a in the region 
m q a < 0.005 implies that there will be little variation 
in lattice simulation of m^a and m p a from our light- 
est quark mass to physical quark mass. Concerned with 
quenched chiral perturbation behavior, we tried chiral 
extrapolations for m v ,m p and from heavier quark 
mass (m q a — 0.0075 ~ 0.05) (these fitting forms make 
sense only in the continuum limit. However, since our lat- 
tice spacing is quite small, we tried various fitting forms 
assuming that there are little modification on the fit- 
ting forms due to the finite lattice spacing effect) and 
compared extrapolated values with our simulated result 
with lighter quark mass (m q a = 0.00125 ~ 0.01). For 
m^ 2 /m q , although a finite lattice volume argument sug- 
gests a singular l/m q behavior, comparison of a fit us- 
ing m w 2 from heavier quark mass with m w 2 from lighter 
quark mass shows that m^/rriq from lighter quark mass 
calculation is less singular than l/m q . In nucleon and p- 
meson, unlike in pion, it is hard to distinguish by chiral 
extrapolation the linear dependence on ?77 T2 , a term ex- 
pected by quenched chiral perturbation theory, because 
1) the statistical fluctuation associated with m^a and 
m p a is larger than that with m^, and 2) the quenched 
chiral perturbation parameter, 6, appears to be smaller 
than that in ref. [] probably due to our larger 6/g 2 . This 
6/<? 2 dependence may be understandable since 6 is re- 
lated to m§/(47r/ 2 ) [] and the asymptotic scaling of m 
can be different from that of /„■. 
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TABLE I. Bare quark mass m, and hadron mass, all 
in lattice units. Under the staggered quark formalism the 
Nambu-Goldstone (NG) pion 7r and non-NG pion 7T2 split at 
0(a 2 ) because of the flavor symmetry breaking, and so do 
the p- mesons p and p2, but the effects are now so small and 
hardly visible. The nucleon is from the even sources which 
give better signal than the corner ones. 

m„a m^a m^ 2 a m„a m„ 2 a mjya 



0.05 0.3845(4) 

0.04 0.3363(4) 

0.03 0.2839(4) 

0.02 0.2266(5) 

0.015 0.1959(6) 

0.01 0.1582(5) 

0.0075 0.1377(6) 

0.005 0.1131(6) 

0.0025 0.0811(7) 

0.00125 0.0576(8) 



0.3890(4) 
0.3394(4) 
0.2868(4) 
0.2284(5) 
0.1962(5) 
0.1577(5) 
0.1394(7) 
0.1121(8) 
0.0850(10) 
0.0612(29) 



0.4196(6) 
0.3767(6) 
0.3322(7) 
0.2882(9) 
0.2661(10) 
0.2434(8) 
0.2347(16) 
0.2229(13) 
0.2137(22) 
0.2122(37) 



0.4198(5) 
0.3770(6) 
0.3316(7) 
0.2878(9) 
0.2664(11) 
0.2417(9) 
0.2367(15) 
0.2225(13) 
0.2122(21) 
0.2117(33) 



0.637(1 
0.568(1 
0.495(1 
0.418(2 
0.380(2 
0.336(1 
0.313(3 
0.293(2 
0.269(3 
0.261(6 




o.oo 



0.25 



0.50 



1.00 



FIG. 1. The nucleon to p mass ratio vs. pion to p mass 
ratio at 6/g 2 = 6.5 for m q a = 0.05, 0.04, 0.03, 0.02, 0.015, 
0.01, 0.0075, 0.005, 0.0025, and 0.00125. Set I is plotted with 
diamonds (o) , and Set II crosses ( x ) . Fit values with only sta- 
tistical errors are shown and no continuum or finite-volume 
correction is made. The lower circle represents the experimen- 
tal value ((mjv/m p , m-Tr/mp) = (1.218,0.182)) and the upper 
circle represents the non-relativistic limit (= (1.5, 1.0)). 
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FIG. 2. The pion effective mass from 12 3 (x), 24 3 (o), 32 3 
(□) and 48 3 (+) wall sources at m q a — 0.01. We see the 
quality of our plateau when all the four agrees. 
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FIG. 3. m%/m q vs. m q a. The curves are fit to a form of 
m%/m q — Co/rn q + C\ + C2in q , to Set I (solid) and Set II 
(dotted). 
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FIG. 4. rriN and m p vs. m^. (x) is nucleon mass from 
Set I and (O) is nucleon mass from Set II. (□) is p mass from 
Set I and (+) is p mass from Set II. The leftmost circles show 
expected tunci and m p a for physical quark mass which is de- 
duced by the nucleon mass and the pion mass from experiment 
and by the current lattice spacing, a. 
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